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Abstract. 

Loeb and Waxman have argued that high energy neutrinos from the decay of pions produced 
in interactions of cosmic rays with interstellar gas in starburst galaxies would be produced with 
a large enough flux to be observable. Here we obtain an upper limit to the diffuse neutrino flux 
from starburst galaxies which is a factor of ~5 lower than the flux which they predict. Compared 
with predicted fluxes from other extragalactic high energy neutrino sources, starburst neutrinos 
with ~ PeV energies would have a flux considerably below that predicted for AGN models. We 
also estimate an upper limit for the diffuse GeV 7-ray flux from starbust galaxies to be O(10 -2 ) 
of the observed 7-ray background, much less than the component from unresolved blazars. 



1. Introduction 

Interactions of cosmic-ray nuclei with interstellar gas nuclei in our galaxy produce 7r°'s which 
decay to produce most of the galactic 7-rays above 0.1 GeV pQ; the decay of the 7r + 's produced 
yields galactic cosmic-ray neutrinos (5). 

The distribution of high energy 7-rays in our Galaxy is related to the distribution of molecular 
clouds and very young hot high-mass stars in OB associations which are short-lived and explode 
into supernovae E] . This association between supernovae which are likely to produce cosmic 
rays and dense regions of molecular gas led to the scenario where interactions between the gas 
and cosmic rays produces the galactic 7-rays via the decay of the ir°. A natural implication 
then would be that starburst galaxies, which are undergoing a phase of extremely active star 
formation would be likely sources of high energy 7-rays [5JE)]- 

Loeb and Waxman (LW) jTj have suggested that such hadronic processes in starbust galaxies 
can produce, in toto, a large enough background of diffuse high energy neutrinos to be observable 
with a very large neutrino detector such as Icecube jSJ. LW then argue that radio observations 
of starburst galaxies imply a lower limit on the cumulative extragalactic neutrino flux from 
starburst galaxies which is within the sensitivity range of Icecube. We derive here an upper limit 
to the cumulative high energy neutrino flux from starburst galaxies which is significantly lower 
than the LW "lower limit" . The diffuse flux of GeV 7-rays from the same processes is found to 
be 0(1O -2 ) of the observed 7-ray background, much less than the component from unresolved 
blazars and more than an order of magnitude below a recent estimate (§]. 



2. Radio Emission and the Neutrino Flux from Starbust Galaxies 

LW start with the observed synchrotron emission from starburst galaxies which is produced 
by relativistic electrons in these sources jlUj . They then make the assumptions that (1) the 
presence of relativistic electrons in these sources implies the presence of relativistic protons, (2) 
the protons lose essentially all of their energy to pion production, and (3), a lower limit to the 
energy loss rate of the protons can then be obtained from the synchrotron radio flux by assuming 
that all of the electrons (and positrons) which are radiating are from pion decay. 

Assumption (1) is a reasonable one which is supported by observations of cosmic rays in 
our own Milky Way galaxy. Assumption (3), viz., the "lower limit" assumption depends on 
assumption (2). However, assumptions (2) and (3) can be questioned because (a) the synchrotron 
radiating electrons may be largely accelerated primaries rather than secondaries related to pion 
production and decay, as is the case in our own Galaxy, and (b) the conditions in starburst 
galaxies are significantly different from those in our own galaxy. In particular, starburst galaxies 
exhibit strong "superwinds" jllj . Such winds have significant dynamical effects and may disrupt 
magnetic fields and drive protons out of these galaxies before they can lose all of their energy by 
interacting with interstellar gas nuclei to produce pions ^2]. In contrast, assumption (2) of LW 
assumes full trapping of relativistic nuclei in the disks of starburst galaxies to the point where 
they only lose energy in hadronic interactions with gas atoms. This is in stark contrast to the 
situation in our own Galaxy (see footnote 1) and is also contradicted by recent observations of 
high redshift galaxies |13j . 

These caveats call into serious question the argument that the radio data can provide a lower 
limit on the cumulative diffuse flux of neutrinos from starburst galaxies. But here we consider 
that assumptions (l)-(3) can be used for obtaining an analytic upper limit for such a flux. We 
also will accept the other estimates which lead to the ratio of injected power of protons to 
electrons at a fixed particle energy, rj p / e ~ 6 and a neutrino luminosity which is then related to 
the local radio luminosity density by 

E%$ V (E V = IGeV) ~ (ct H /An)t[Af(dL f /dV)] f=lAGHz (1) 

where tu is the age of the universe and £ = 3 is an evolution factor which takes account of 
the fact that starburst galaxies were more numerous in the past 0. LW take the local energy 
production rate per unit volume at a frequency f = 1.4 GHz to be ~ 10 28 ' 5 W Mpc -3 . Let us 
reexamine this value for f(dLf/dV)]f=i,4GHz- 

The local 1.4 GHz energy production rate has been derived by LW by making use of the 
correlation bewteen GHz and far infrared (FIR) emission in galaxies given by Yun, Reddy and 
Condon (YRC) ^U]. YRC use the data on IRAS galaxies to derive the local infrared luminosity 
density at 60 /mi to be 2.6 x 10 7 L Q Mpc" 3 . This total power density is then used by LW to 
obtain the 1.4 GHz power density via the relation given by YRC. The key difference between the 
result derived here and that obtained by LW is in chosing how to interpret the paper of YRC. 
YRC state that less than 10% of the local FIR luminosity density is contributed by luminous 
IR galaxies with Lfir > 1O 11 L ; this is the component which includes the starburst galaxies. 
(Figure 11 of YRC yields an estimate of ~ 6%.) We therefore take the local contribution from 
starburst galaxies alone at 60/^m to be < 2.6 x 10 6 L Q Mpc~ 3 . Consequently, the component of 
the local radio luminosity density from starburst galaxies is <&iaghz < 10 27,5 W Mpc~ 31 

There is a higher relative fraction of the energy input from the higher relative number of 
starburst galaxies at higher redshifts. The fraction of the FIR background, k(Az) contributed 
by galaxies in different redshift ranges, Az, is obtained from Ref. [Tl]. We multiply k(Az) by the 

1 One might ask "Why not consider "normal" galaxies with lower FIR luminosities and add them in to estimate 
a higher neutrino flux?" However, in normal galaxies like ours, cosmic rays lose only a small fraction of their 
energy via hadronic interactions, contrary to assumption (2). Also galactic PeV cosmic rays have a much steeper 
spectrum than the E~ 2 assumed by LW. 



fraction of the FIR background contributed by starburst galaxies in different redshift ranges, 
£(Az), to estimate the mean fraction of the total FIR background contributed by starburst 
galaxies. Estimates for k and £ are shown in Table 1. Using the results from Table 1, we 
estimate that 23% of the observed FIR background integrated over redshift is from starburst 
galaxies. 



Table 1: Relative contributions to the v Starburst Galaxy Flux (see text). 



Redshift Range (Az) 


k(Az) 14 




Reference for £ 


to 0.2 


10% 


< 10% 




0.2 to 1.2 


68% 


~ 13% 


pu 


>1.2 


22% 


~ 60% 


m 



3. Observability of High Energy Neutrinos from Starburst Galaxies 

The upper limit on the radio flux from starburst galaxies obtained above can be used to obtain 
an upper limit on the neutrino flux from starburst galaxies by using equation (1) of LW. One 
then finds that the neutrino background energy flux from starburst galaxies < 2 x IO -8 GeV 
cm -2 s -1 sr Such a flux would be undetectable above the atmospheric background neutrino 
flux, even if equation (1) is assumed to be valid when extrapolated to 300 TeV and even granting 
all of the assumptions made by LW. 2 



Table 2: Neutrino Energy Fluxes (GeV cm 2 s 1 sr 1 ) 
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Table 2 shows a comparison of the upper limit on the flux from starburst galaxies given here 
with the atmospheric neutrino flux and with approximate model predictions of neutrino fluxes 
7-ray bursts (GRB) and active galactic nuclei (AGN) along with detector array sensitivities. 
It can be seen from this table that at 100 TeV none of the extragalactic sources proposed will 
dominate over the atmospheric foreground. The present upper limit on the diffuse neutrino 
energy flux below 1 PeV from AMANDA-II is ~ 8.8 x IO" 8 GeV cm -2 s" 1 sr" 1 in the 10 TeV to 

1 PeV energy range [18] . The full Icecube detector array is expected to push down to a sensitivity 
of ~ IO -9 GeV cm" 2 s -1 sr -1 in the energy range 100 TeV < E u < 100 PeV after several years of 
observation. Under the extreme assumption that the primary cosmic ray spectra in all starburst 
galaxies are as hard as E~ 2 up to energies O(10 PeV), PeV neutrinos from starburst galaxies 

2 Even if we make a second extreme assumption that 100% of the IR galaxies at redshifts greater than 1.2 are 
starburst galaxies, we would still predict a neutrino flux < 3 x 10 -8 GeV cm -2 s _1 sr -1 . If we use the fast 
evolution model of Ref. I25| . we would obtain a similar upper limit. 



may be detectable just above the projected sensitivity of Icecube. However, as can be seen from 
Table 2, above 1 PeV the AGN models predict fluxes which will be significantly larger than the 
the atmospheric foreground (expected to be < 3 x 10~ 9 ) GeV cm -2 s~ x sr _1 at 1 PeV), as well 
as the fluxes from 7-ray bursts (GRB) and starburst galaxies. 

4. Observability of Diffuse 7-rays from Starburst Galaxies 

In a follow-up paper to LW, Thompson et al. [5] have estimated the contribution of 7r°-decay 
7-rays from starburst galaxies to the observed 7-ray background in the GeV energy range. Using 
an E~ 2 primary spectrum they get estimates of E 2 &(E) fluxes for both 7-rays and neutrinos of 
3 x 1(T 7 GeV cm' 2 s" 1 sr" 1 . 

Our estimated starburst galaxy 7-ray flux for the same E~ 2 primary spectrum assumed in 
Ref. [H] is ~ 2 x 1CT 8 GeV cm -2 s -1 sr _1 for 7-ray energies less than ~ 10 GeV. This is 0(1CT 2 ) 
of the observed flux of ~ 1.4 x 1CT 6 GeV cm~ 2 s^ 1 sr _1 determined by the EGRET group [2S] 
and is therefore unobservable. Above ~ 10 GeV the background spectrum will steepen owing 
to absorption from pair production interactions with the extragalactic ultraviolet background 
radiation 24 . Stecker and Salamon have shown that the bulk of the observed background can 
be produced by unresolved blazars [23]. Note, almost all of the observed extragalactic GeV 
7-ray sources are blazars; no starburst galaxies have been observed. 
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